Method to estimate the parameters of a propagation channel 

R A rfCdU OT TND Q gTTTF INVENTION 
1. Field of the Invention 

The invention relates to a method for estimating the impulse response of 
a propagation channel, especially its parameters such as the delay values, the 
directions of arrival as well as the values of the complex amplitudes associated with 
these parameters, with an a priori knowledge of the signal. 

It is applicable, for example, to the estimation of the parameters of a 
finite impulse response filter which can be written, though not necessarily so, in 
specular form, namely a filter that can be written in the form of weighted D.racs. 

The invention can be used to estimate propagation channels in the field 
of radiocommunications but also, generally, it can be applied to any signal filtered by 

a finite impulse response linear filter. 

In a transmission system, especially one using radio waves, a transmitter 
sends out a signal in a transmission channel to a receiver. The signal that is sent 
undergoes amplitude and phase fluctuations in the transmission channel. The signal 
received by the receiver consists of copies of the transmitted signal that are 
temporally shifted and modified. The fluctuations of the signal and the shifts give 
rise to a phenomenon known to those skilled in the art as inter-symbol interference. 
The interference arises especially out of the law of modulation used for transmission 
and also from the multi-path propagation in the channel. 

The received signal generally arises out of a large number of reflections 
in the channel. The different paths taken by the sent signal thus lead to different 
delays in the receiver. The impulse response of the channel represents the totality of 
the fluctuations to which the sent signal is subjected. 

The estimation of the propagation channel in a radiocommunications 
system is useful in several respects, some of which are indicated here below as 
examples. 

. The demodulators generally require knowledge of the channels ,n order to 
remedy the harmful effects that they have caused, 

• a second point of interest is that of urban or extra-urban localization, for example 
the principle of emergency localization using the number "91 1 » in the United States 

• Finally, the knowledge of the propagation channels can also serve for the use of 
smart antennas at reception as well as transmission. 

2. nescri ption the Prior Art 
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There are various techniques known in the prior art for estimating 

propagation channels and their parameters. 

For example the document by R. Rick and L. Milsteil, "Performance 
acquisition in mobile ds-c4ma systems", IEEE Trans on Communications, Vol: 45 
5 (No: 1 l):pp: 1466-1476, November 1997, proposes a search for propagation delays 
by using a bank of non-coherent detectors. The results are proposed for multi-path 
channels in the presence of Doppler phenomena and inter-cell and intra-cell 
interference. 

The document by R. Rick et L. Milsteil, "Optimal decision strategies for 
10 acquisition of spread-spectrum signals in frequency selective fading channels" in 
IEEE Trans, on Communications, Vol: 46 (No: 5): pp: 686-694, May 1998, discloses 
an optimal decision rule based on the outputs of the correlators proposed in the 
document referred to here above. A single-user technique of this kind is limited by 
interference in the case of multiple-users. 
15 There also exist known ways of using rectangular shaping filters, for 

example by the method described in one of the following documents: 
. E. Strom, S. Parkvall, S. Miller, and B. Ottersen, "Propagation delay estimation 
in asynchronous direct-sequence code-division multiple access systems", IEEE 
Trans on Communications, Vol: 44:pp: 84-93, January 1996 
20 • S Parkvall, "Near-Far Resistant DS-CDMA Systems : Parameter estimation and 
Data Detection", PhD thesis, Royal Institute of Technology Stockholm, Sweden, 

1996. , . . 

. S. E. Bensley and B. Aazhang, "Maximum likelihood estimation of a single user s 

delay for code division multiple access communication systems", Conf. Information 

26 Sciences and Systems, 1994. 

In the case of shaping filters with a duration greater than a chip time, 

these different methods are no longer suitable. 

Algorithms for the combined estimation of angles of arrival and of the 
differential delay times, on known and received signals, based on sub-space 
30 techniques have been proposed, for example in the document by P. Gounon, 
-Analyse spatio-temporeile haute resolution a 1'aide d'une antenne active", (High 
Resolution Space-time Analysis Using an Active Antenna ) Traitement du S.gnal 
(Signal Processing), Vol. 1 1 (No. 5), pp. 351-360, 1994. 

The document by A. J Van der Veen, M. C. Vanderveen, et A. J. Paulraj, 
35 "Joint angle and delay estimation using shift-invariance properties", IEEE Sig. Proc 
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Letters, Vol.4 (No.5): pp. 142-145, 1997. discloses methods for the estimation of the 
physical parameters of propagation by means of methods based on sub-spaces. 

However, such methods suffer from a deterioration of performance 
characteristics once the impulse responses of the propagation channels are correlated. 
5 This situation occurs especially when the complex amplitudes do not vary with 
sufficient speed on the covariance matrix of the impulse responses estimated in terms 
of the least-error squares by means of the signal transmitted. 

Figure 1 shows the different techniques of maximum likelihood. 

A maximum likelihood method has been proposed, for example, in one 

10 of the following references: 

• J. Grouffaud, "Identification spatio-temporeJIe de canaux de 
propagation a trajets multiples", (Space-time Identification of Multi-Path 
Propagation Channels), PhD thesis, Ecole Normale Superieure de Cachan, June 
1997. 

• M. Wax and A. Leshem, "Joint estimation of delays and directions ot 
arrival of multiple reflections of a known .signal.", IEEE Trans, on Signal 
Processing, Vol: 45(No: 10):pp: 2477-248, October 1997. 
but it does not deal with the MIMO (Multiple Input Multiple Output) context. 

The document by P. Graffouliere, "Methodes actives spatio-temporclles 
large bande" (Active Wideband Space-Time Methods), published in « Techniques et 
performances. Applications En Sonar » (Techniques and Performance. Sonar 
Applications), PhD thesis, INPG. March 1997, also discloses a method of estimation 
based on maximum likelihood but the studies on performance deal only with the case 
of a single source or of several distinctly separate sources. A similar .study is 
disclosed in the document N. Bertaux , "Contribution a Tutilisation des methodes du 
Maximum de Vraisemblance en traitement radar actif ' (Contribution to the Use of 
Maximum Likelihood Methods in Active Radar Processing), PhD thesis, Ecole 
Normale Superieure de Cachan, January 2000, lor active radar applications in the 

case of single sources. 

30 SUMMARY Q F THF. INVENTION 

The present invention relates to a method that can be used especially to 
estimate the parameters of the propagation channel by working on the correlated 
signals, in selecting a certain number of samples and in searching for the values of 
the delay parameters and/or directions of arrival, for example, which would enable 

35 the most efficient reconstruction of the signal received. 
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The invention also relates to a method that integrates pulse compression 

techniques such as pre-processing. 

The invention relates to a method for the estimation of one or more 
parameters of a propagation channel with a priori knowledge of the signal in a 
system comprising one or more sensors 



The method of the invention comprises the following steps: 
rrelating the signal or signals x(t) received by the sensor or sensors with 

a known signalctl)?*^^,^^ 

• sampling said^iglfel-aiie^correlation at a sampling period Te and selecting 

a number of samples per concatenation>>»^^ 

. determining at leas tone parameter^hepropagation channel such as X 
and/or 8 which enables the most efficient reconstmc^^th^signals received by 
using a maximum likelihood method. 



According to one embodiment, the characteristics^ the system of 
sensors are known and comprise for example: 
. a correlation step with a known signal c(t) equal to 1 , 

. the signals received on the antenna being expressed in the form X=S(6,T)h+B, 
and 

. the estimates of the parameters as t and 9 being expressed in the following form: 

0 t T = Mgmin||n£(0,T)X|| a 
= arg min 

where n^s is the projector orthogonal to the image generated by the 

column vectors of S(t, 0). 

The characteristics of the system of sensors arc for example the response 

to the antenna. 

It may comprise a step for determining the complex amplitudes h of the 
impulse response of the propagation channel from the estimates of the estimated 

parameters as x and 6. 

According to another embodiment, the characteristics of the system ol 

sensors are not known, and the method comprises for example: 

• a step for the correlation of the signals received by the network of sensors 
with a known signal c(t) equal to I , 



. a step in which the signals received are expressed in a concatenated form 
Y = y(t)a+N where y(T) is equal to the convoluted product of the unit matrix In 
and the matrix S^Kt^Ar *(r#/)J, and a contains the responses of the 

paths of the different users, 
5 • a step for the estimation of the delay vector t from 

r = argmin|ini(r)Y|f 
= argm'n/WY'niMY) 

where 17^ is the projector orthogonal to the image generated by the line vectors of 

to y(x). 

The step of correlation of the signals is performed with a signal c(t) 
different from 1. The characteristics of the system of sensors are known and this 
correlation step comprises a step for the estimation of the parameters T and 6 from: 
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0 t - argminX w R^ , ni(0,T)X lJ 



where n ± ^j^^t) {&{0,T)K?*{0,Tj)- 1 ^{B,r)^ 1 
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The invention also relates to a device for estimating one or more 
parameters of a propagation channel with a priori knowledge of the signal in a 
system comprising one or more sensors. 

The device comprises at least: 
• a device adapted to the correlation of the signal received by the sensor or 

25 sensors s(t) with a known signal c(t), 

. a device adapted to the selection of a number of samples of the signal 

obtained after the correlation step, and 

- a device adapted to the determining of the parameters of the channel by a 

maximum likelihood method. 

The method can be applied, for example, in applications of the MIMO 

(Multiple Tnput Multiple output) or SIMO (Single Input Single Output) type. 

The invention in particular has the following advantages: 
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. in the cooperative MIMO context, the estimator, within the framework of a 
known antenna, makes use of a parametric model of the signals received on the 
antenna, unlike the deterministic maximum likelihood method which estimates each 
sample of the signal without any a priori assumptions, 

. in limiting the field of analysis in certain cases of application, it leads to a 
reduction of the size of the data processed, and thus results in faster processing. 

RRTEP DESCRIPTION OF TP * DRAWINGS 

Other features of the invention shall appear from the following 
description which provides examples of implementation given by way of an 
illustration, with reference to the appended drawings, of which: 

. Figure 1 shows different approaches of estimation of parameters using the 

Maximum likelihood method 

• Figure 2 gives a schematic view of a first alternative embodiment of the 
method according to the invention in which the characteristics of the reception 

antennas are known. 

. Figure 3 is a drawing of a second alternative embodiment in which the 

characteristics of the reception antennas are not known, 

. Figure 4 gives a schematic view of another variant of the method 
comprising a preliminary step that can be applied to the methods described in figures 

2 and 3 here above. 

morf. DETAJLE D nFsrRIFTlON 

In order to provide a clearer understanding of the object of invention, the 
description is given by way of an illustration that in no way restricts the scope of the 
invention in the field of mobile radio communications, in a cooperative context m 
which the receiver has certain information available to it, for example the 
knowledge, to a certain extent, of the signal sent (a priori knowledge ot the signal) 

The receiver is constituted for example by N identical and 
omnidirectional sensors. The signals s(t) received by each of the sensors are 
regularly sampled at the rate T e . This example does not prevent the application ot the 
method of the invention to other types of receivers. 

Before giving a detailed description of the different alternative torms ot 
implementation of the method, we shall recall a few points concerning the modelling 
of the signals and shall state a few assumptions enabling the implementation ot the 
method. 

Modelling of signals received on the reception system 
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When there are U users (sources), the U signals s u (t) with 1 * u S U are 
propagated through different channels before being received by a set of sensors 
indexed k. The impulse response of the channel between the user u and the sensor k 
is referenced h u k (t ? T) where t is the instant of arrival of the signal and t is its delay. 

If we call h u k (t, t) the vector with the dimension (Afxl) whose 

ii. i ku t* f\ vuith I < k < K and K the number of sensors, the 
components are the values h k(t, wun i _ *■ - 

models of the received signals can be written in the form ( 1): 

x(t)=Xs„(T) *h U (t,t)+b(t) where 



b(t) is the noise vector with the dimension ( N x 1 ). 
1 0 Specular Assumption 

Time-related specular assumption 

The specular model is a widely used model of channels for which it is 
considered that the paths are perceived by the antenna as being temporally discrete 
paths. 

The transfer function of the channel therefore becomes (2): 



hU(t,T> £ h>p u (t))S(t-T p u(t)) 

with 

• P u the number of propagation paths of the user u , 

h U (T p "(f)) the (//xl) vector containing the product of the responses ot the 

antenna to the path p multiplied by the complex amplitude of the path. 
The received signal can be written in the case of the specular channel (3) 

x(o-tf:w ))A ^ <(0)+b(0 

Calibrated antenna 

If the antenna is calibrated, assuming that the signals reaching the 

25 antenna have a negligible angular dispersion and the assumption that the signals are 
narrowband signals is met, the received signal is modelled by (4): 



11-1 p-1 

where 

# 9 u is the direction of arrival of the path p of the user u, 
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h U (T o "(0) is the complex gain of the path p of the user u received on the first 



sensor. 



The previously made assumption (concerning the specularity of the 
direction of arrival) is not restrictive in the case of a base station in height. If the 
5 specularity is not sufficient, the model may be modified by adding a lobe w.dth m 
the direction of arrival which will be estimated by the method proposed ,n the 
invention by adding the parameters of this lobe to the parameters to be estimated. 

The different variants of the method according to the mvention use the 
maximum likelihood method, for which various approaches to the estimation of 

10 parameters have been given in figure 1 . 

First alternative embodiment: calibrated antenna 

A first alternative embodiment of the method according to the invention 
uses a maximum likelihood method adapted to the signals sent, when these signals 
sent are partly known, and the antenna response for all the directions of arrival ,s 

15 known. " . 

Figure 2 shows a block diagram where the processing block I, 

comprising for example a microprocessor, adapted to implementing the method 
according to the invention, receives firstly the different signals received by the 
antennas 2i as well as a known signal corresponding to the signal sent. 

At output of the processing block, the method according to the mvention 
has at least estimates of the direction of arrival 9 and delay vectors x which, for 
example, are used in order to deduce therefrom the impulse response of the channel 

h(t) for a given user or for all the users. 

The signals received on the antenna array are delayed, expanded and, 
possibly frequency-transposed versions of the signals sent by the users present in a 
cell (or propagation channel). These signals are therefore expressed as a linear 
combination of a set of known basic functions whose parameters are to be extracted. 

The method is presented in the general context of parametric estimation 
of the angles 9 and of the delays x for any antenna whatsoever constituted by N 
identical and omnidirectional sensors. The signals arc sampled regularly at the rate 
Te It is assumed that the number of paths are known and that, if this is the case, the 
number is determined by using a method of detection of the number of sources by 
which is possible to obtain the number of paths, described for example in the 
following references. 



• J. J. Fuchs, "Estimation of the number of signals in the presence of unknown 
correlated sensor noise", IEEE Trans on Signal Processing, Vol 40 (N° : 5) pages 
351-360, 1994; 

. J.Gouffraud, P.Lanzabal and H.Clergeot, "Methodes a Haute resolution pour le 
traitement d'antenne et l'analyse spectrale" (High-resolution methods for antenna 
processing and spectral analysis) in « Selection d'ordre de modele », Editions 
Hermes, February 1998, Chapter 6; 

. H. Akaike, "A new look at the statistical model identification", IEEE Trans. On 
Automatic Control, Vol 19 (No. 6), December 1974, 

• M Wax and T. Kailath, "Detection of signals by information theoric criteria" de 
M.Wax et T.Kailath, IEEE Trans.on ASSP, Vol 33(N°2), pages 387-392, April 

1985. 

The steps implemented in the algorithm according to the invention 
remain valid when the number of sources is overestimated. 
Modelling 

The signal x(t) (equation (4) received by the array of sensors 2i is 
sampled at the rate T e during a period t=[T e K e T e ] with Ne being the number of 
samples or number of observations during which the complex amplitudes ot the 
impulse response for a user u and a path p. h u p <T« p (t)) are considered to be 

constant. . 

These N e observations with the dimension Nx\ are concatenated tn the 

vector X with the dimension N e Nx\: X =[ x(T 9 ) T ,-,X(N e T e ) ] (5) where the 

exponent T designates the transpose. 

We define the vector B with the dimension N e Nx\ contaming the 
concatenation of the samples of the noise during the same period by (6): 

B=[b(T 6 ) T ,-,b(N e T e ) T l T 

Finally, s w (x p u ) is taken to denote the vector with the dimension N e x l 
corresponding to tne signal sent by the user u, delayed by (V and sampled: 
s"((V)=£ •s u (T e -(Tp U ),-,s u (NeT e -(V')] (7) 

The vector of the concatenations of the signals received can therefore be 
written as follows (8): 

u~l p=t 
where ® is Kronecker's product 



Let S(9,x) be the matrix N e N*P ( P=±p" ) of the delay-direction vectors 

u=i 



defined as follows (9) : 



where: 

- 0 is the P-sized real vector of the directions of arrival . 

- i is the P-sized real vector of the delays, (P being the number of 
paths) 

The equation of the signals received on the antenna can therefore be 
written as: X=S(9,t)h + B (10) where h is the Pxl complex vector of the amplitudes. 

This formulation has the value of expressing the observations as a hnear 
model that is a function of the complex amplitudes h and of the noise. 

It must be clearly noted that this is a MIMO (multiple input multiple 
output) modelling, grouping together all the signals of all the users. 

It is not necessarv to know all the signals sent. For example, m mobile 
communications, sequences of known signals are sent periodically, resulting in a 
known signal. In the radar or in radio altimetry, for example, the transmitted pulse ,s 

known. 

Log-likelihood 

The complex amplitudes h are considered to be unknown but 
deterministic. In order to simplify the explanation of the invention without, however, 
limiting it, the signals s u (t) forming S(G.t) are considered to be perfectly known, 
apart from the parameters fct This leads to an observation model (10) where only 
the additive noise vector is random, with a Gaussian probability density, represented 

by (11): 

F{X|*h>.r> = - s;5 L^exp(-i(X-S (e ,r)h)'(X- S (< > .r,h,) 
The log-likelihood of the observations is therefore (12) : 



L (X|** h. 0, r) - W) - £ IX - S(fl, r>h||' 
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Estimator of the parameters C\ h, 8, x 

Through the separability of the parameters a\ h, 6, x, the expression of the 
log-likelihood (12) authorizes the analytical determining of the vector of the complex 
amplitudes h and of the power of the noise o 2 as a function of Q and X. These 
estimates are both obtained by cancelling the derivatives of the log-likelihood 
respectively with respect to a 2 and with respect to h. Their expressions are 
respectively (13) and (14): 



10 and 



h = (S f (9, r)S(e, r))~ l &(&, r)X 



By substituting the values of h and of a 2 in the log-likelihood (12), the 
1 5 expression gets simplified and the estimates of 9 and x are given by ( 1 5) 

0, r = arg min X - S(*. r) (*<•. t)S(9, r» " S» (*, r)X 



In letting n* denote the projector orthogonal to the image generated by 
the column vectors of S(8, t) (16) 

nJ^rj.I-W.rXBttftrW.rM-^.r) 



The estimates are the solutions of (17) 

min||n^(0,T)X|| 2 



0,9 = argmi 



These estimated values are then used respectively in the expressions 
(13) and (14) to determine the power of the noise and the impulse response of the 
propagation channel. 



The criterion introduced in (17) is referenced hereinafter in the text by: 

L( 0,r)=»X'Il£«? f T)X 

Implementation of the method according to the invention in the case of a non- 

calibrated antenna 

Figure 3 gives the schematic view of an exemplary implementation of the 
method in the case where the characteristics of the receiver or the antenna are not 
known. The antenna is formed by several sensors, whose response is not known, and 
is often designated by the expression "unstructured antennas". 

The antennas 2i receive the signals that are transmitted to a processing 
block 14 adapted to implementing the steps of the method, this block itself rece.v.ng 

a known signal. . 

The processing block, which is equipped tor example writ a 

microprocessor programmed accordingly, delivers at least the estimated value of the 
delay vector X and the estimated value of the impulse response of the propagation 
channel. 

Modelling 

The received signal is expressed in the form of the equat.on (4) and ,s 
sampled at the rateTe on a period of time t = ITe, NeTe] in which the complex 
amplitudes of the impulse response of the propagation channel h p (tp (0) may be 
considered to be constant or substantially constant. These Nc samples w.th a 
dimension N*l are concatenated in a vector Y with a dimension NeNx. such us 
(18): 



Y = 



xi(r,),...,^7;r) ,.- ^(r,),...,^(Mr t ) 



Where xi(nTe) is the sample n of the sensor i. 

Similarly, the NeNxl sized vector N contains the concatenation of the no.se samples 



(19). - 



N - 



h(TA....MNJrA. MT e )....^(N.T,) 

. \ J s v 



Where bi(nTe) is the noise sample n of the sensor with an index i. 
The signal Y can be written (20) 
Y = \i/(T)a+N 



where the matrix (N^xNP) y is built from the different values s"(x p u ) as follows 
(21): 

= I N ®S(X) 

where S( T > is the matrix (W„xP) such as (22): 

a contains the responses of the sensors to the different paths of the different 
users (23) : 



a = 



with: 



The complex amplitudes are always assumed to be constant during the phase 

of estimation of the parameters. 

The modelling of the received signals given by the equation (20) provides us 
with an affine model . As above, this expression takes account of the signals of all 

the users. 

Log-likelihood 

The complex amplitudes a are considered in this case as being unknown but 
deterministic . The signals that come into play in 4>(t) are assumed to be known apart 
from the variables T . This leads to a model where only the noise N is random, with 
Gaussian components and with a probability density (25) : 

P {YkV.r} - js^exp (-> «Y - *(r)a,.(Y - *<r)a») 
The log-likelihood is then given by (26) 

i( Y|.',a, T ) = - h l|Y " * (rQ)f 

Estimator of the parameters o*, a, x 

As in the case of the maximum likelihood adapted to the signals received 
with a known antenna as explained here above, the analytical expressions of the 



A 
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complex amplitudes a as well as the power of the noise* may be expressed as a 
function of the set of parameters of the delay vectors t. 

This can be obtained for example by respectively cancelling the deri vauves of 
the log-likelihood with respect respectively to a and to a . 

We obtain the expression (27) 

and (28) 

S = (»»(r)*(r))- i '••MY = * # M Y 

By replacing a and J by their estimates, the log-likelihood gets simplified 
and the estimator of t is (29) 

r = argmin||Y - *<r) (••(rWr))-' *'(r)Y| 

If n«t(t) denotes' the projector that is orthogonal to the images generated by 
the line vectors of \|/(x) (30): 

ni(r) _i-*M(*'W*M)- *'M 

The estimate of T is given by (3 1 ) 

f = org mm ||niM Y ll 

- argntintr(Y'niM Y ) 

The two alternative embodiments expired here above deal with signals 
S ampled at output of the sensor or again for a signal corrected with a known signal 
havine a value 1 . The matrices to be processed have a s,ze NeNxP. 

According to another alternative implementation of the method which can be 
applied^ the c JL of a calibrated antenna as we., as that o, a non-calibrated antenna, 
the method according to the invention uses correlated signals. 

Figure 4 shows a block diagram in which a correlator ,0 receives, firstly, the 
sianals received by the antennas and. secondly, the known signal. 

8 The correlator .0 has a number of outputs identica. to the number o, mpu s 
each signal sift) that reaches an antenna being collated by a known signal c(t>. The 
d £2 corre ated signa.s s.(„-ciC. where i is the index of the centre a. , en 
Emitted to a device 1 1 adapted to the selection of a number o. samples, espcca.ly 
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in order to work on a time interval defined beforehand- The signals that have been 
correlated and belong to the chosen time interval are then transmitted to a processing 
unit 12 comprising, for example, a microprocessor suited to estimating the values 9 
and T, when the characteristics of the antenna are known or the value of i for an 
5 antenna having unknown characteristics. From these values, it is possible to 
determine the estimates of the impulse response of the propagation channel. 

The correlator 10 may be of the analog or digital type, the sampling steps 
being performed before or after the correlation step. 

The processing block also receives the sent signal, itself correlated for 

10 example by means of a correlator 13. 

When a known time interval is used, it is useful to know the time of reception 

at the block as well as the length of the channel. 

The signals considered in the method correspond to the signals output after 
correlation by the filter matched with the transmitted signals. 

Method for an antenna whose characteristics are known. 
In this part, the same assumptions are used as in the previous part explained 
in the case of an antenna for which the parameters are known. We shall again 
consider the signals to be known and shall examine only the case of the user CO in 
order to simplify the presentation. This model can immediately be extended to the 
20 MIMO (Multiple input Multiple Output) context. 

Modelling 

By correlating the signals received at each sensor with the signal of the user 
CO, we get (32) : 

= fEf;a(W(r;)^-r;) + b(0j*^'(-0 

- 22 «(W ( 5< " ( * - r " tt) * sW( " 0) + (b(0) * 

By limiting the data processed to the part "containing" the propagation 
channel at output of the correlation, the complexity is greatly reduced because the 

size of the matrices used is diminished. 

By concatenating in the filtered samples at output of the sensors on a time 
interval Nl, the signal X can be modelled as follows (33) : 



35 



10 



15 



20 



X = *(0,r)h+b 



where 



1? 
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and 



with 



and *,(t,e) = W-T;)*s"(-t))*m 

The noise b 40 follows a normal law of covariance OrRb where the matrix 
R h =Rc m Iw • R/ is the covariance matrix of the vector (34): 

The likelihood of the signal X" is given by (35): 



1 V(X"-*(*.T)h-) , R^ , (X--*^.'-)h«) 

F(X*|*(*.r),W-,4,Rb) = ^ <y ^|R b | €XP ' T 



whence the log-likelihood (36): 

L(X w |<&(0, T), h w , <rf , Rb) - -^' iV T/ "t 

By cancelling the derivatives of the log-likelihood with respect to o* and with 
respect to h"\ we obtain the estimators of a* and h w (37). 

n 



*2 — 



and 



R -i (X w _ $( 0, r)h w ) 
TViN 



By replacing a* and h<° by their expression in (36) we obtain the estimator sought for 
the space-time parameters 8 and T (38); (39) 
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where 



0,t = argminX-R^^i^.rJX 4 - 
ni -!-*(<?, r) (*'(*. fl-)Hi l Wr))" , » t («.r)Rj l 



The estimated space-time parameters are used for example to determine the 
35 estimate of the impulse response of the propagation channel. 



Non-calibrated antenna 

In this part, we shall look at a reduced maximum likelihood method with an 
unknown antenna. The approach used will be identical to the preceding one. We 
shall give only the model of the signal, the log-likelihood and the estimators. 

Modelling 

The received signal is filtered by j^C-f) (40) : 

xf(r) - ftf h> W |f(t "* ,+b( T^ H 

By concatenating the samples at output of the sensors and then the vectors 
obtained, we get (41): 



where 



and S(r) = [Si , • • • , Sp-J 

With S p = [S p (T e ) l --- f S P (N l T t )\ T 

Sp(0 - 5 W (« - r p u ) * ^(-f) 



and 



The noise n u follows a standard law of covariance afR, where the matrix 

The likelihood of the signal Y* 0 is given by (42): 

POr-IS(r), of, o>, Rn) = ^^77 jgj ~P V 7 
whence the log-likelihood (43) : 

L(Y»]S(t), a w , 4, R„) = -Wf/V log(Traf) - log(|R„|) - i I R-n ' (Y" - S(t)cO 



By cancelling the derivative of the log-likelihood with respect to a\ and with 
respect to a w , we obtain the estimators of a\ (43) and of a m (44): 



-=2 _ 

a T — 



Rn * (Y W - S(t)o?) 



^ a"« (S(r)m I S(T-)) - l S(r)ltf X" 

By replacing in (42) 6$ and <*" by their expression, we find an 
estimator for the parameter T (45): 

9 = axgmin VR^nfti-JV 

T 

where n£ =1- S(r) (S(r)B^ ! S(r))" 1 S(r)R; 1 

The estimated value of the delay t is then used for example to estimate 
the value of the impulse response of the propagation channel. 

The different alternative modes of implementation of the method 
according to the invention, can be applied in many fields of signal processing. In 
communications systems, it is often necessary to know the propagation channel in 
order to demodulate the signals transmitted. 

In view of spectral congestion and the growing demand from users, 
future systems will have several users communicating simultaneously in the same 
frequency band. Consequently, the techniques proposed here in the MIMO context 
have a direct application. 

Similarly, the technique of filter estimation with reduction of complexity 
can be applied in the MIMO and SIMO (Single Input Multiple Output) context. 

These techniques are also applicable to radar processing. They can be 
used in a multiple-target or single-target context with echos 

Similarly, in radio altimetry, several echoes of the sent signal may be 
received and all of them have to be estimated. 
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Without departing from the framework of the invention, the invention 
may take account of the Doppler effect. For this purpose, one method consists for 
example in estimating the speeds by minimizing a multidimensional function that 
integrates all the parameters of the model, other than the complex amplitudes and the 
noise power. 

The method according to the invention can be applied especially in a 
UMTS telecommunications system comprising an array of antennas. 



